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ABSTRACT 

An analysis of an experiment to determine the precision, the effects 
of varying amounts of environmental data, and the feasibility of the use 
of historical environmental data at the NAVTORPSTA, Keyport, Washington, 
Dabob Bay facility is presented. The analysis verifies and extends the 
results of a preliminary investigation of the experimental data. The ex- 
periment consisted of tracking a 75-kHz source at six depths at each of 
six horizontal ranges. The data were analyzed utilizing an isogradient 
ray- tracing program. The results showed that the precision of the range 
position as determined by the ranging system is dependent only on the 
rise- time of the 75kHz signal; that the use of less environmental data 
and historical data as the environmental input for position calculations 
have been shown to cause only a minor degradation in the accuracy of the 
ranging system. 



4 



TABLE OF CONTENTS 



I. INTRODUCTION 8 

II. DABOB BAY TRACKING FACILITY 9 

A. DABOB BAY FACILITY 9 

B. RANGE DESCRIPTION 10 

III. EXPERIMENTAL PROCEDURE 12 

IV. RANGE DETERMINATION TECHNIQUES 14 

V. PRELIMINARY ANALYSIS 38 

A. ANALYSIS 18 

B. RESULTS 18 

C. CONCLUSIONS 19 

VI. PRECISION OF DATA 20 

VII. ENVIRONMENTAL EFFECTS ON RANGE PRECISION 22 

VIII. HISTORICAL ENVIRONMENTAL ANALYSIS 26 

IX. CONCLUSIONS 29 

X. COMMENTS AND RECOMMENDATIONS 31 

FIGURES 32 

TABLES 54 

APPENDIX A. ISOGRAD 64 

APPENDIX B. ISOVEL 67 

APPENDIX C. ANGCAL 71 

APPENDIX D. STUTRACK I 76 

BIBLIOGRAPHY 85 

INITIAL DISTRIBUTION LIST 86 



5 



LIST OF ILLUSTRATIONS 



1. DABOB BAY RANGE 32 

2. EXPERIMENTAL AREA 33 

3. HYDROPHONE ARRAY 34 

4. TILT CORRECTION 35 

5. RAY PATH REFRACTION 36 

6. PHI versus THETA, by Station 37 

7. PHI versus THETA, by Depth 38 

8. COMPARISON OF SOUND VELOCITY PROFILES, as derived by the 

Equal-Change-in-Velocity Method 39 

9. COMPARISON OF SOUND VELOCITY PROFILES, as derived by the 

Equal-Change-in-Gradient Method 40 

10. COMPARISON OF SOUND VELOCITY PROFILES, as derived by the 

"Eyeball" Method 41 

11. DABOB BAY SOUND VELOCITY PROFILE, TYPE 1 42 

12. DABOB BAY SOUND VELOCITY PROFILE, TYPE 2 43 

13. DABOB BAY SOUND VELOCITY PROFILE, TYPE 3 44 

14. DABOB BAY SOUND VELOCITY PROFILE, TYPE 4 45 

15. DABOB BAY SOUND VELOCITY PROFILE, TYPE 5 46 

16. DABOB BAY SOUND VELOCITY PROFILE, TYPE 6 47 

17. DABOB BAY SOUND VELOCITY PROFILE, TYPE 7 48 

18. DABOB BAY SOUND VELOCITY PROFILE, TYPE 8 49 

19. DABOB BAY SOUND VELOCITY PROFILE, TYPE 9 50 

20. DABOB BAY SOUND VELOCITY PROFILE, TYPE 10 51 

21. DABOB BAY SOUND VELOCITY PROFILE, TYPE 11 52 

22. DABOB BAY SOUND VELOCITY PROFILE, TYPE 12 53 



6 



LIST OF TABLES 



I. STANDARD ERROR OF ESTIMATE, THETA 54 

II. STANDARD ERROR OF ESTIMATE, PHI 55 

III. STANDARD ERROR OF ESTIMATE, R 56 

IV. STANDARD ERROR OF ESTIMATE, caused by Environmental 

Description by Velocimeter 57 

V. STANDARD ERROR OF ESTIMATE, caused by Environmental 

Description by the Equal-Change-in-Velocity Method 58 

VI. STANDARD ERROR OF ESTIMATE, caused by Environmental 

Description by the Equal-Change-in-Gradient Method 59 

VII. STANDARD ERROR OF ESTIMATE, caused by Environmental 

Description by the "Eyeball" Method 60 

VIII. STANDARD ERRORS OF ESTIMATE OF APPROXIMATED PROFILES, 
as derived from Historical Profile 4 between the 

Surface and 200 feet 61 

IX. STANDARD ERRORS OF ESTIMATE OF APPROXIMATED PROFILES, 
as derived from Historical Profile 4 between 100 

and 300 feet 62 

X. COMPARISON OF STANDARD ERRORS, caused by the Approximation 

of Sound Velocity Profiles by the Equal-Change-in-Velocity 
Method 63 



7 



I. INTRODUCTION 



An acoustic experiment was conducted to provide a basis for the 
comparison of various ray-tracing techniques which can be utilized to 
determine the three-dimensional position of a source with respect to a 
fixed receiver. The experiment was conducted at the Dabob Bay facility 
of the Naval Torpedo Station (NAVTORPSTA) , Keyport, Washington on 21 
March 1974 by members of the Ray Trace Committee of the Range Study 
Group at the Naval Postgraduate School, Monterey, California. 

The purpose of the experiment was to determine what factors, both 
instrumental and environmental, influence range accuracy and precision. 

A preliminary investigation of two of the fifty-seven sets of data 
collected [Ref. 1] led to these tentative conclusions: 

(1) The precision of the position location is essentially independent 
of the range from the array and depends only on the rise-time of the 
75-kHz signal. 

(2) There is no significant difference in the accuracy of results ob- 
tained by either isovelocity or isogradient ray-tracing techniques. 

(3) The isogradient method appeared to need less environmental data 
to obtain a given precision. 

This thesis was undertaken to verify and to extend the results of 
this preliminary investigation by testing the isogradient method of 
ray-tracing at all stations and depths studied and by determining the 
amount of environmental data required to maintain accuracy and by ascer- 
taining if historical data rather than real time data may be utilized. 
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II. DABOB BAY TRACKING RANGE 



A. DABOB BAY FACILITY 

A facility for the three-dimensional tracking of surface ships, 
submarines, and torpedoes is located at the Dabob Bay facility of the 
Naval Torpedo Station, Keyport, Washington. Chosen because of its 
favorable oceanographic features and its proximity to the NAVTORPSTA, 
Keyport, Dabob Bay has been used for torpedo testing since 1949. It 
consists of a 75-kHz primary tracking system enveloping nearly the 
entire bay providing for underwater tracking in a volume approximately 
30,000 feet by 4,500 feet by 600 feet deep. 

Dabob Bay, a branch of the Hood Canal, is in a deep depression ad- 
jacent to the Olympic Mountains. The Quilcene River flows into the 
northernmost part and the Dosewallips River into the southernmost. The 
sides of the bay are precipitous and predominately rocky, while the 
bottom is mud. Annual precipitation averages about 15 inches and the 
average snowfall is about 15 inches. Winds generally blow along the 
length of Dabob Bay in a North-South direction. Normal wind velocities 
are 5 to 15 miles per hour. Tidal levels in Dabob Bay range from about 
-5 to +15 feet referenced to mean low-low water (MLLW) . 

Water temperatures change markedly from season to season and even 
from day to day. Salinity varies considerably with depth as fresher 
surface water from terrestrial runoff overlies more saline water. The 
variation in salinity is also quite seasonal, reaching a maximum in the 
spring. The main contributor to sound velocity changes is water tem- 
perature variation and the resulting sound velocity profile is 
representative of the temperature gradient. 
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B. RANGE DESCRIPTION 

The basic underwater tracking system is comprised of three compo- 
nents: hydrophone array, computer system, and vehicular instrumentation. 

The shipboard transducer emits a 75-kHz acoustic pulse in synchronism 
with a master clock at the computer site. This pulse is detected by 
each of the four hydrophones in the array and mixed with a local oscil- 
lator having a frequency unique to that hydrophone. 

The outputs of these oscillators are preamplified and fed to a 
multiplexer, which sums the signals and transmits the composite signal 
via underwater cable to the computer site, where it is processed to de- 
termine the in-water transit time to each hydrophone (Tx, Ty , Tz, Tc) . 

1. Hydrophone Array Configuration 

The hydrophones in the NAVTORPSTA Range are arranged in a 
short-base line system in groups of four, each one located on one of 
four adjacent corners in an imaginary 30-foot cube, thereby defining 
the orthogonal coordinate system in which the measurements were made. 
These hydrophone arrays use buoyant spheres which exert an upward force 
in excess of two thousand pounds to keep them upright. The array is 
not free to rotate about the Z-axis , but may tilt as much as two degrees 
from the horizontal in strong currents. This tilt is measured by servo- 
pendulum transducers whose specified resolution is 0.00833 degrees, 
accuracy 0.025 degrees, and linearity 0.05 per cent over the full scale. 

The range consists of six arrays (Fig. 1) aligned along a line 
450 feet east and parallel to the range axis, spaced at 6000-foot in- 
tervals, with the 00 array being the farthest north. The range axis is 
191° 18* 14" T. 
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2 . Computer System 



The computer system, a Scientific Data Systems (SDS) Model 920, 
is installed on Zelatched Point and consists of three main subsystems: 
the signal-processing subsystem, the data-collection subsystem, and the 
computer subsystem. The signal-processing subsystem, the link between 
the hydrophones and the data-collection subsystem, processes the multi- 
plexed signal received from the arrays, discriminates against unwanted 
signals, and determines the in-water transit time to each hydrophone. 

The data-collection subsystem establishes the master clock timing by 
which all timing is computed. It also prepares all the array signals 
for the computer subsystem, which calculates the position of the tracked 
object, prints the tracking data, and records the data on magnetic tape. 
For the real-time examination, the data were plotted on X-Y plotters. 
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III. EXPERIMENTAL PROCEDURE 



A 92.5-dB (re 1 //bar) acoustic source with approximately ten pounds 
of ballast was suspended from the stern of a sound boat to depths of 25, 
50, 75, 100, 150, and 200 feet at six horizontal ranges (stations) from 
900 feet to 4700 feet from the 02 hydrophone (Fig. 2). The first sta- 
tion was located approximately 900 feet north of the array, and 250 feet 
east of the range axis. Successive stations were located parallel to 
the range axis; the second being directly over the array, the third 600 
feet past it, and then at about 1500-foot intervals to about 4500 feet, 
beyond which signal- to-noise problems made further examination impossi- 
ble. The source emitted one 1.3 millisecond 75 kHz pulse per second for 
sixty to ninety seconds at each depth at each station. 

The horizontal position of the boat on the range was found from the 
range 1 s Auto tape system, a portable commercial microwave system used 
exclusively for ship tracking that yields accuracies of + 3 feet. The 
position of the acoustic source with respect to the Autotape antenna was 
determined from the length of cable let out, the horizontal distance 
from the antenna to the stern of the boat, and the ship’s heading. 

Tidal data were recorded and used for correcting depth coordinates to 
MLLW. The array is known to be 585 feet below MLLW, so that the depth 
of the source below the surface and the tide measurement allows calcula- 
tion of the vertical distance between the source and the array. All 
measurements were made with respect to the array’s coordinate system, 
although the array’s orientation and location on the range were known, 
making it possible to transform this data to a true position on the 
range. The X-axis of the 02 array was 25° 24’ 36.0” T. 
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A NAVTORPSTA digital velocimeter that has an accuracy to +1.0 ft/sec 
and resolution of + 0 . 7 ft/sec was used to measure the sound-velocity 
profile at each station. NAVTORPSTA instrumentation at Dabob Bay meas- 
ured and recorded the acoustic travel times from the source to each of 
the four hydrophones in the array to six decimal places for every source 
depth at each station. The three-dimensional position of the source 
with respect to the array was calculated by the NAVTORPSTA NUTRAK III 
computer program. 

The data were gathered on 21 March 1974 during a neap tidal period 
and while the bay surface was smooth and glassy, although there had been 
ripples earlier in the day. The day was bright and sunny with scattered, 
thin, high overcast, the temperature was 55°F, and no noticeable wind was 
present. The shiny stainless steel transducer disappeared from view when 
about eight feet deep, indicative of thick biological concentrations in 
the water. Rudimentary drift measurements of boat movement relative to 
the surface water were taken by throwing paper into the water and watch- 
ing its relative motion with respect to the boat, which indicated a value 
of approximately 0.3 ft/sec. 



IV. RANGE DETERMINATION TECHNIQUES 



A. TECHNIQUES 

The three-dimensional position of an acoustic source with respect to 
a receiver may be readily determined by computer analysis. Two ray- 
tracing techniques which are easily adaptable for high-speed digital 
computation are: (1) an isovelocity technique in which the water mass is 

divided into layers of isovelocity water; and, (2) an isogradient tech- 
nique in which the water is divided into isogradient layers. Both 
techniques utilize Snell T s law, which describes the refraction of sound 
rays in a medium of variable velocity. 

The isovelocity method is currently utilized in the NAVTOR.PSTA pro- 
cedure NUTRAK III. STUTRACK I was developed in the preliminary 
experimental analysis. It consists of both an isovelocity and an iso- 
gradient method to be utilized as comparators with NUTRAK III. Both 
programs require as inputs the velocity profile, the transit time and 
angle of arrival of the acoustic wave at the geometrical center of the 
array. 

B. PREPARATION OF DATA 

The required inputs for the computer analysis must be derived from 
the information available at the array, which is the transit times of the 
acoustic pulse from the tracked source to the four hydrophones of the 
receiving array. The geometry involved in the calculation of the raw 
array coordinates is shown in Figure 3. The expressions for the raw 
array coordinates are: 
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RAWX = (T c + T x ) (T c - T x ) (1) 

2d 

2 

RAWY = -§r- (T + T ) (T - T ) (2) 

2d c y c . y 

and 

2 

RAWZ = (T + T ) (T - T ) (3) 

zd c z c z 

where 

c = the average sound velocity, ft/sec; 
d = the distance between hydrophones, 30 ft; 

T^, Ty, T^, = travel times to respective hydrophones, sec. 

A correction factor must be applied to correct for the array receiver 
integration time (pulse width discrimination), 800 /is . The relations now 
become: 

o 

RAWX = — — (T' + T' - K) (T' - T') (4) 

zd c x c x 

c 2 

RAWY = -Z— (T' + T' - K) (T ’ - T') (5) 

zd c y c y 

and 

2 

RAWZ = -2— (T' + T' - K) (T' - T') (6) 

zd c z c z 

where 

-3 

K = twice the integration time, 1.6 x 10 s. 

T 1 , T 1 , T 1 , T ! = the measured transit time, sec. 
x y z c 

These raw coordinates are referenced to a non-horizontal reference 
plane due to the array’s freedom to tilt in the horizontal. The array is 
instrumented to sense the angles the array axis makes with the horizontal, 
XTILT and YTILT. 

The correction of the RAW coordinates is a matter of geometry, as 
shown in Figure 4. The corrected coordinates become: 



CORX = RAWX - RAWZ [SIN(XTILT)] 

CORY = RAWY - RAWZ [SIN(YTILT)] 

CORZ = RAWZ -i- RAWX [SIN(XTILT)] + RAWY [ SIN (YTILT) ] 



(7) 



( 8 ) 

(9) 

In a homogeneous bodv of water, the minimum time path between an output 
transducer and a receiving hydrophone is a straight line. When the speed 
of propagation is spatially varying due to temperature and salinity vari- 
ations, refraction occurs and the minimum time path is no longer a straight 
line (Figure 5) . Corrections must therefore be made to determine the true 
position of the source. 

From Snell f s law, a minimum time path is computed, working back in 
time from the array center in the direction of the "apparent" position and 
tracing the acoustic ray through each velocity layer until the measured 
time is consumed. The vertical angle of entry at the array, computed from 
the tilt-corrected array coordinates is 



A = 
o 



SIN 



CORZ 2 

CORX 2 + CORY 2 + CORZ 2 _ 



1/2 



( 10 ) 



This entry angle is assumed to be equal to the entry angle for the 
minimum time path. This assumption is good if the source is far enough 
so that the wave fronts are planar over the dimensions of the array, and 
the speed of sound is constant over the dimensions of the array. 

The travel time, T, to the origin is computed from 



T 



T 

c 




(ID 



where 



S = t (CORX) 2 + (CORY) 2 + (CORZ) 2 ] 1/2 



( 12 ) 



and 
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( 13 ) 



222 1/2 
R = [ (CORX + 15) + (CORY + 15) + (CORZ + 15) ] ' 

c 

where R^ is the slant range to the array center hydrophone, and S is the 
slant range to the array center. 

The entry angle, A^, the travel time, T, along with velocimeter data 
describing the sound velocity profile may now be utilized as inputs to 
the computer program. The outputs will provide the three-dimensional 
position of the source with respect to the receiver. 
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V. PRELIMINARY ANALYSIS 



A. ANALYSIS 

The analysis of the experimental data was to have determined the 
accuracy of NUTRAK III as compared with both modes of STUTRACK I. The 
Autotape system was to form the reference for this determination, but its 
precision was much less than expected. The analysis was, therefore, pri- 
marily devoted to the determination of the precision of the various 
methods of ray- tracing. Stations 1-25 and 6-50 were chosen for prelimi- 
nary investigation since station 1 was closest to the array and station 
6 was farthest. The effects of varying the amount of environmental data 
was also investigated. 

B. RESULTS 

NUTRAK III showed good point to point continuity at station 1 for all 
coordinates, while the coordinates at station 6 were very erratic. A 
transformation of the coordinates to a spherical coordinate system showed 
that the variations of phi versus theta to be approximately the same at 
both stations. 

Both modes of STUTRACK I were compared to NUTRAK III and in both 
cases no significant differences were noted when the most structured 
sound velocity profile was utilized. As the number of linear segments 
defining the sound velocity profile was decreased, the isovelocity tech- 
nique provided consistently poorer results. Using the isogradient 
technique to compute the source’s position resulted in only a slight 
degradation as the number of segments in the profile was decreased and was 
considered adequate in precision. 
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C. CONCLUSIONS 



Comparison of the two stations investigated showed station 6-50 to 
have greater spatial variability. The transformation to spherical co- 
ordinates shows the difference is not present in the angular display as 
would be expected if the time errors were the important source of vari- 
ability. Improved range accuracy can only be obtained by improving the 
accuracy of the time measurements or increasing the baseline of the 
hydrophone array. The isogradient method of STUTRACK I provides greater 
precision than the two isovelocity techniques when less environmental 
data was utilized. 
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VI. PRECISION OF DATA 



The data collected at Dabob Bay were tested at all stations to verify 
that the results of the preliminary analysis were generally applicable. 

If the precision was independent of range and depth the standard devia- 
tion of the source coordinates as calculated from the transit times at 
the array should be zero with a small standard deviation about the mean. 
Factors which could cause variation are the rise-time of the 75-kHz source 
(as discussed in the preliminary analysis) , and the actual drift of the 
source transducer caused by the currents which exist in the experimental 
area. 

Analysis of the precision of the ranging system must remove the ef- 
fects of drift from the calculations. Regression analysis on the 
components of the calculated horizontal position of the source provided 
an equation for the displacement of the transducer as a function of time. 

A second-order regression equation was considered to adequately describe 
this displacement as it accounts for the velocity and acceleration com- 
ponents of the drift. 

Each set of data was processed to determine the calculated components 
of the source position. Each horizontal component was fitted to a simple 
second-order curvilinear regression equation with time being the indepen- 
dent valuable and the corresponding component as the dependent variable. 
The value of the component evaluated from the regression equation for each 
time was taken as the estimated value at that time. By taking the dif- 
ference between the estimated and the corresponding calculated values the 
effects of drift are removed from the computations. The mean depth was 
taken to be the simple arithmetic mean. The standard deviation, standard 
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error of estimation, of these differences between the estimated and cal- 
culated coordinates is a measure of the precision of the ranging system. 
Tables I-II1 show the standard errors for the theta, phi, and radius 
components of the source position for all stations and depths. Figures 
6 and 7 are scatter diagrams of the standard errors of phi versus theta 
by station in Figure 6 and by depth in Figure 7. Both Figures show that 
the standard errors with the drift effects removed are essentially random 
with no dependence on station or depth. 
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VII. ENVIRONMENTAL EFFECTS ON RANGE PRECISION 



The effect on range position accuracy and precision due to varying 
the amount of environmental information was investigated by varying the 
number of linear segments which were used to describe the sound velocity 
profile. Decreasing the number of segments which describe a particular 
profile while maintaining the required precision would result in faster 
processing time for the ranging data. 

Three methods of decreasing the number of segments in a profile were 
studied: (1) equal changes in gradient; (2) equal changes in velocity; 

and (3) n eyeball tr fitting. All three methods were compared to determine 
which would suffer the least degradation in accuracy and precision as 
the profile was simplified. In the experiment, environmental data was 
provided by a velocimeter which described the existing sound velocity 
profile in 64 to 67 discrete points in 600 feet. Any profile described 
adequately with fewer points would significantly reduce data processing 
time. The equal- change-in-gradient and velocity methods utilized in this 
investigation were developed by the author for this analysis. 

The equal-change-in-gradient method determines a set of linear seg- 
ments by locating the depth and velocity points which, if connected by 
straight lines, would cause a given percentage change in the gradient. 

The first three points from the bottom of the most detailed sound velo- 
city profile are fitted with a straight line by the method of least 
squares. The slope of this line is used as the initial reference gradient. 
The gradient of the linear segment in the original profile between the 
first point used to determine the reference gradient and the next point 
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in the profile is computed. This gradient is compared with a test value 
given by the reference gradient plus a designated change from 5 to 30%. 

When the computed gradient exceeds the test value, the prior point in 
the profile is considered as the end point in that segment. The points 
thereby included in this segment are fitted by the method of least squares. 
The end point of the segment which has been processed is now the first 
point used in the determination of a new reference gradient as described 
above and the next linear segment is determined by repeating the process. 
The velocity at the highest depth of the initial segment does not equal 
the velocity at the lowest point of the second segment after processing; 
therefore, an adjustment must be made to cause the two segments to inter- 
cept at this mutual depth. This adjustment consists of making the sound 
velocity of both profiles the mid-point of the calculated velocities at 
that depth. This process is continued until the surface is reached. 
Appendix A is a listing of the computed program which utilizes this 
method for decreasing the number of segments which describes a sound 
velocity profile. 

The equal-change-in-velocity method utilizes a constant percentage 
change in velocity between the minimum and maximum velocities in the pro- 
file and the minimum and maximum velocities to decrease the number of 
segments in the profile. The minimum and maximum sound velocities of the 
original profile are located; a fixed percentage (5 to 30%) of the dif- 
ference between them is used to determine the velocity at each point in 
the new profile. The surface sound velocity is incremented by the given 
change in velocity and the original profile is searched until the depth 
of the incremented velocity is located or the minimum or maximum velo- 
cities occur. The method determines the direction of the change in 
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velocity from the original profile. Once the new point in the profile is 
determined, its velocity is incremented and the sequence, as described 
above, continues until the bottom of the original profile is reached. 

The original sound velocity profile is now described by a lesser number 
of linear segments. Appendix B is a listing of the computer program 
which utilizes this method. 

The ,, eyeball ,, method is simply the author’s estimate of the best- 
fitted segments which describe a sound velocity profile. It is derived 
by estimating the segments which provide the best description of the 
original profile in the minimum number of segments. Then by decreasing 
the number of possible segments, estimates are made of the best-fitted 
profile with the decreased number of segments. 

Thirteen curves describing the sound velocity profiles at both sta- 
tions 1 and 6 were determined using the above methods, six by the equal- 
change-in-velocity method, three by the equal-change-in-gradient method, 
and four by the "eyeball" method. Figures 8-10 show the change in the 
characteristics of the profile at station 1 due to its description by 
each of these methods. These profiles were used to study the effect on 
source position due to varying the amount of environmental data. The 
position of the source at each depth for stations 1 and 6 was calculated 
for each profile. Tables IV-VII show the resultant standard errors of 
estimate from these calculations. At both stations the standard error 
of the theta and phi components remained constant as the amount of envi- 
ronmental data was varied. The standard error of the slant range changed 
with each profile. Excluding the results for the "eyeball" method at 
station 6, the average standard error for all methods was 0.181 feet at 
station 1 and 0.317 feet at station 6. The average standard errors by 
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method for stations 1 and 6 were: (1) equal-change-in-velocity , 0.160 

and 0.343 feet; (2) equal-change-in-gradient, 0.220 and 0,318 feet; and 
(3) "eyeball", 0.172 and 1.586 feet. The equal-change-in-velocity method 
appears to cause the least degradation in positional precision with the 
least environmental data at these two stations. 
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VIII. HISTORICAL ENVIRONMENTAL ANALYSIS 



The possibility of using historical environmental information rather 
than real-time information as the sound velocity input for computing 
position has as its major advantage the lessening of the cost and time 
required prior to performing any range operation. Historical information 
along with its reduction by the methods discussed in Chapter VII would 
decrease total processing time and, thereby, increase the efficiency of 
the ranging system. This hypothesis is only valid if the use of this in- 
formation does not degrade the accuracy of the positional calculations. 

The objective of this section, therefore, is to analyze the feasibility 
of utilizing historical environmental information in computing position 
with precision as the measure of effectiveness. 

Figures 11 thru 22 show the twelve sound velocity profiles which de- 
scribe the environmental conditions which exist at Dabob Bay. The 
profiles were based on three years (1971-1973) of environmental informa- 
tion collected at the facility [Ref. 2]. Each profile type shows the 
minimum, maximum, and median (mid) profiles which are characteristic of 
the particular type. These historical sound velocity profiles were used 
in this analysis. 

In order to perform the analysis, the inputs for the range positioning 
algorithm, STUTRACK I, had to be determined based on the historical pro- 
files. A simulation was developed to generate the transit time of an 
acoustic pulse from a source to the array, the angle of entry of the 
acoustic wave with respect to the center of the array, and the apparent 
horizontal distance of the source from the array. The simulation 
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generates these parameters by searching for the entry angle in which a 
ray would enter the array from a known position of a source. Once the 
entry angle is located, the transit time for that angle is determined. 
This simulation, as with STUTRACK I, utilizes Snell’s law as its basis 
and performs the same calculations as STUTRACK I, but in reverse order. 
Appendix C is the computer listing of the subroutine, ANGCAL, developed 
for the simulation. 

The simulation computes the input data for STUTRACK I at a fixed 
horizontal position from the array. This position was chosen to provide 
a horizontal range from the array similar to that of station 6 in the 
experiment. Fifty vertical points at four-foot intervals from the sur- 
face to 200 feet were used as the vertical component of the source 
position. The use of this scenario was chosen to provide data for the 
analysis under the conditions that were found to cause the least preci- 
sion of source position in the experiment, station 6. 

Each of the historical profiles was first described in 60 ten-foot 
layers from the surface to 600 feet, and this description was used as the 
environmental input for the simulation. The transit time, entry angle, 
and apparent horizontal range for the 50 positions were used to compute 
the three-dimensional position of the source with STUTRACK I. Each 
historical profile was then reduced using the equal-change-in-velocity 
method at five to thirty percent change intervals, providing approxi- 
mations to the original profile. Each of these simplified profiles was 
used as the environmental input to STUTRACK I and the source position 
was calculated with the transit times, entry angles, and horizontal 
ranges from the original profile. 
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The position of the source as calculated from the reduced profiles 
was compared to that of the original profile. The standard error of es- 
timate for each of the reduced profiles was determined for the 36 
profiles describing the twelve historical profiles at Dabob Bay. A 
comparison of the minimum and maximum profiles of a particular type with 
the mid profile was also performed. This comparison was designed to show 
the effects of using a single profile to describe the environment when 
the actual profile existing at that time is within the range of a parti- 
cular sound velocity type. Table VIII shows the results of these 
comparisons for Profile Type Pour. All profile types showed similar 
results . 

In conducting the above comparisons, it was noticed that much larger 
errors occurred at shallower depths. Each profile type shows its great- 
est variation near the surface due to water runoff, solar insolation, 
and seasonal temperature variation. Therefore, it was decided to perform 
the same analysis as above, but with the data points below the point 
where large variations in the gradients of the minimum and maximum pro- 
files of a sound velocity type occurred. The calculations were performed 
between 100 and 300 feet. Table IX shows the results of these calcula- 
tions for Profile Type Four as before. The values, as calculated from 
data points between 100 and 300 feet, generally, showed less standard 
error than the values calculated using the same number of points between 
the surface and 200 feet. Table X shows this comparison. 
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IX, CONCLUSIONS 



A. PRECISION OF DATA 

The analysis of the experiment showed conclusively that although the 
magnitude of the spatial variation grows as the source’s range from the 
array increases, the variation in the angular description of source 
position is independent of range. This is consistent with the results 
of the preliminary investigation. It is concluded that the increase in 
the spatial variation of the range position discrimination of the 2.5-MHz 
clock in measuring the arrival time of the 75-kHz source signal used in 
the experiment. This uncertainty can cause a maximum error of 3.3 ^sec 
in the arrival time which is one-fourth of the period of the 75-kHz sig- 
nal, causing an error of approximately one foot at station 1 and 4.5 
feet at station 6. 

B. ENVIRONMENTAL EFFECTS ON RANGE PRECISION 

The investigation of precision of range position due to the amount of 
environmental information used in the calculations has shown that the 
method of describing the existing profile is not critical. Of the three 
methods of reducing the number of segments describing the existing sound 
velocity profile, the equal-change-in-velocity method appears to produce 
results with more precision than either the equal-change-in-gradient or 
'’eyeball* 1 fit methods. Describing a profile in terms of a 20 percent 
equal change in velocity between segments, usually six or seven segments, 
produces the best precision with the minimum amount of environmental 
information required for calculations. This result is independent of 
depth and horizontal range. 
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C. HISTORICAL ENVIRONMENTAL ANALYSIS 



The analysis of the twelve profile types which describe the acoustic 
environment at Dabob Bay show that the use of historical information 
rather than real-time environmental data is indeed feasible. The use 
of a profile which has been reduced by a 20 percent change in velocity 
in the equal-change-in-velocity method provides the most precision for 
all velocity profile types studied. The use of the mid profile to de- 
scribe any profile falling within the range of a particular type causes 
very little degradation in precision. Therefore, it is concluded that 
the use of historical data is feasible in the determination of the source 
position and only a small number of segments are required to describe the 
existing sound velocity profile. 
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X. COMMENTS AND RECOMMENDATIONS 



This thesis has verified the preliminary conclusions reached in 
Reference 1, that the precision of the range position in calculations is 
a function of the rise-time of the 75-klIz acoustic signal used in the 
experiment. The angle of entry of the signal to the array is thus the 
critical factor in range precision. Increasing the accuracy of time 
measurement and decreasing the error due to rise-time (increasing the 
array baseline) will improve precision. 

The historical profiles may be segmented as described in this thesis 
and catalogued. The sound velocity profile which actually exists on the 
range could be sampled daily with a velocimeter. This profile can then 
be fitted to an historical profile type which describes it and the ap- 
propriate catalogued sound velocity profile used as the environmental 
input for any range positioning calculations to be performed that day. 

As noted in the preliminary investigation, the sum of the source 
depth and the individual calculated Z component at each time differed 
from the known array depth. The range of this difference was +4.8 to 
-18.8 feet with the average being 3.1 feet less than the array depth. 
Verification of the recorded tidal levels with the predicted, tides at 
the time the experiment was conducted showed no discrepancies of this 
magnitude. 

Further study should be conducted to determine the accuracy of the 
ranging system at Dabob Bay and Nanoose and to provide a greater con- 
fidence level in the use of less environmental data in range position 
calculations. The apparent inaccuracy in the Z component should 
receive special attention. 
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RAUZ SIN(XTILT) 



The SIN(XTILT) is positive when the X transducer is 
above the horizontal plane passing through the C 
transducer, and negative when the X transducer is 
below the horizontal plane passing through C. 

FIGURE H. TILT CORRECTION 
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FIGURE 5 . RAY PATH REFRACTION 
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TANDARO ERRORS OF ESTIMATE OF APPROXIMATED PROFILES 
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STANDARD ERRORS OF ESTIMATE OF APPROXIMATED PROFILES 
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STANDARD ERRORS CALCULATED BETWEEN 100 AND 300 FEET 



COMPARISON OF STANDARD ERRORS 
CAUSED BY THE 

APPROXIMATION OF SOUND VELOCITY PROFILES 
BY THE EQUAL CHANGE IN VELOCITY METHOD 



STANDARD ERROR 
RANGE OF CALCULATIONS 



PER 

CENT 

CHANGE 



0 - 200 FT 
PHI ! R 
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(FT) ,i 
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0.03l! 


1 . 684! 
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i.9aa! 


0.031*. 
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15 


• 0©016* 


3.306! 


O.Oli! 


1.45l! 




20 
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5.17?! 


0.037! 


1 .681 ! 




25 


« 0 • Q^<L • 


6.086! 


0.034! 


1.961*. 




30 


• 0 • 0 3 • 


6.213! 


0.012*. 


2 . 048 ! 


TABLE X 
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EGUAL CHANGE IN GRADIENT METHOD 
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EQUAL CHANGE IN VELOCITY METHOD 
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SUBROUTINE ANGCAL 
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STATION NUMBER, NUMBER OF DEPTHS TO BE STUDIED, 
NUMBER OF POINTS IN THE SVP, AND CLOCK HCJR DATA WAS 
TAKENtFORMAT - A2, 3110). 
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